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Certain analysis methods have emerged recently that claim to reveal 
"higher harmonic flows" in more-central A-A collisions at the RHIC and 
LHC. But pQCD calculations describe the same structures quantitatively. 

PACS numbers: 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 

1. Introduction 

I review methods for obtaining "higher harmonic flows" from correlation 
data in the context of a two-component model of pt spectra and correlations 
and pQCD calculations of spectrum hard components. The main focus is 
on the mechanism for the same-side (SS) 2D peak in angular correlations. 

2. Two-component model of p-p collisions 

To understand the structure of high-energy A-A collisions we must first 
develop an accurate p-p reference. The two-component model of measured 
200 GeV p-p spectra and correlations includes a soft component - projec- 
tile nucleon dissociation (proton fragmentation) - and a hard component - 
large-angle-scattered parton fragmentation dominated by 3 GeV gluons. 

Figure [T] (first) shows spectrum hard components from 200 GeV p-p col- 
lisions on transverse rapidity yt = \n[(pt + rrit) / m.,^] which provides balanced 
visual access to low- and high-pt structure [Ij . We observe that soft and hard 
components scale with Ug and respectively (n<j is the soft component of 
multiplicity rich) and are thus easy to distinguish down to 0.5 GeV/c. 

Figured] (second) shows a followup study of two-particle correlations on 
{ytiUt) showing distinct soft and hard components. Figure [1] (third, fourth) 
show corresponding soft and hard components of angular correlations. The 
data are minimum-bias, with no trigger condition. Structure can be split 
into same-side (SS) and away-side (AS) on azimuth, with like-sign (LS) and 
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Fig. 1. First: ?/t spectrum hard component vs 7ic/j,, Second: Correlations on (yt, j/t), 
Third: Correlations soft component, Fourth: Hard component (jet structure). 



unlike-sign (US) charge combinations. The combined systematics of spec- 
trum and correlation hard components compel a (di)jet interpretation [21 13]. 

3. Correlation and spectrum systematics for Au-Au collisions 

Given a well-defined p-p reference what happens to the two-component 
model in A- A collisions? Extensive study of Au-Au spectra and correlations 
at 62 and 200 GeV establishes centrality and energy dependence [H [5] . A- A 
centrality is measured by mean participant path length v = 2Nbin/Npart 
obtained from the Glauber model, u is used only as a geometry parameter 
(the N-N cross section is maintained at the 200 GeV value for all cases). 

Figure [2] (first, second) shows angular correlations for most-peripheral 
and most-central Au-Au collisions. The data extend down to 95% central ~ 
N-N collisions. The most-peripheral data correspond well with p-p results. 
Data from the peripheral region are required to establish a Glauber linear 
superposition (GLS) reference corresponding to transparent A-A collisions. 




Fig. 2. First: 200 GeV Au-Au correlations for 85-95% central. Second: 0-5% 
central. Third: Same-side 2D peak amplitude. Fourth: Away-side ID peak ampli- 
tude. 



Figure [2] (third, fourth) show results from 2D model fits that extract 
all information from the data [5|- A simple parametrization includes three 
principal model elements: (a) SS 2D peak, (b) AS ID peak on azimuth and 
(c) nonjet (NJ) azimuth quadrupole strongly evident in more-central Au-Au 
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collisions - a "third component." The remaining structure (also modeled) 
is conversion electron pairs, Bose Einstein correlations and a ID peak on r/A 
(soft component) that decreases to zero amplitude by mid centrality. The 
plots show the SS 2D peak and AS ID peak amplitudes for two energies. 

The minimum-bias SS 2D peak is monolithic, with no separate "ridge" 
feature. The AS ID peak on azimuth representing all aspects of transverse 
momentum conservation is dominated by parton scattering to back-to-back 
jets. The centrality dependence reveals a sharp transition in certain jet- 
related structure properties, a large change within one 10% centrality bin in 
the slopes of SS and AS peak amplitudes (relative to GLS references) and 
the SS 2D peak r/ width. The deviations from GLS for the SS 2D and AS ID 
peaks are quantitative. Both remain consistent with jet-related structure. 

A two-component study of Au-Au pt spectra for identified hadrons was 
also conducted [1] . The soft component maintains a fixed form (participant 
dissociation remains unchanged). Although the hard component changes 
substantially with centrality it is still described quantitatively by pQCD j6]. 

4. SS 2D peak and Fourier series and NJ quadrupole 

We now address the title issue: how to find or not find "higher harmonic 
flows" in A-A correlation data. In the past two years it has become popular 
to project all 2D angular correlations onto ID azimuth, fit the projection 
with a Fourier series and interpret all series terms as "harmonic flows." We 
now relate such procedures to the two- (or three-)component model. 

Figure [3] (first) shows the Fourier coefficients for a Gaussian with r.m.s. 
width £70^ given by Fm{(T^^) = y/JfH a^^(i-x.^{-m'^a1^/2) iSj- The 
hatched bands labeled SS and AS correspond to the typical azimuth widths 
of SS and AS jet-related peaks. The AS ID peak is described by a dipole 
(m = 1). The SS 2D peak requires several multipoles (Fourier series terms). 




Fig. 3. First: Fourier coefficients vs width for ID Gaussian, Second: SS 2D peak 
for 0-5% central Au-Au , Third: SS peak ID projection with Fourier components. 

We next consider the azimuth structure of 2D correlations for 0-5% 
central Au-Au collisions. For that centrality Aq{2D} ~ [9j. The 2D data 
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histogram minus the fitted AS dipole leaves only the SS 2D peak, as shown 
in the second panel. The fit residuals are consistent with statistics. The 
only source of higher harmonics in those data is the SS 2D (jet) peak. 

Figure [3] (third panel) shows the ID projection of the SS 2D peak 
(bold histogram) and its Fourier components (bold dotted curves) with 
amplitudes described by 2Ax{SS}{b) = (6)]G[c7,,^ (6), Ar/]^2D(^) = 

2pQ{b)v'^{S S} (b) . Factor G represents projection of the SS 2D peak onto 
ID azimuth [8j. X represents various multipoles (e.g., dipole D, quadrupole 
Q, sextupole S, octupole O) with 2m poles for cylindrical multipoles, and 
SS denotes multipoles derived from the projected SS 2D peak. 

Figured] (first panel) shows quadrupole amplitudes for various V2 "meth- 
ods." The NJ quadrupole Aq{2D} is derived from 2D model fits with sim- 
ple systematics described by AQ{2D}{b) = 0.00A5Nf,inR{^/sNN)e'^^ (solid 
and dashed curves) [9]. Aq{2D} is related to the Vm by AQ{2D}{b) = 
po{b)v2{2D}{b) with po{b) = nc/i/27rAry dn(.h/2iTdr]. The NJ quadrupole 
is by definition distinct from SS and AS jet-related 2D correlation structure. 




Fig. 4. First: Quadrupole amplitudes for various V2 methods. Second: Quadrupole 
amplitude vs collision energy, Third: Jet-related vs nonjet quadrupoles compared. 

Figure S] (second) shows energy variation proportional to log(-^/sjviv) 
factors R, R': nucleon collectivity below 13 GeV and a novel (QCD?) phe- 
nomenon above 13 GeV carried by a small fraction of the final state [9l 110). 

Figure H] (third) shows the relation among jet-related structure, the 
NJ quadrupole and published V2 data. The basic relation is Aq{2}(6) = 
AQ{2D}{b) + AQ{SS}{b). AQ{SS}ib) is derived entirely from fitted SS 2D 
peak properties, including strong width variations with centrality. We find 
that jet-related and nonjet quadrupole terms from a 2D model fit sum to 
published V2{2} data from the two-particle cumulant "method" [111. The 
dotted curves in the first and third panels are the same, combining the 
measured NJ quadrupole [9] and SS 2D jet peak [5] trends. The {2} and 
{EP} methods are statistically equivalent [12]. Note the small error bars 
within the open squares for V2{2} data in the first panel. The published 
uncertainties have been multiplied by 5 to make them visible in the figure. 
Deviations from Aq{2D} data (upper solid curve) are tens of error bars. 
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5. Azimuth multipoles vs pQCD jets 

We now consider examples of finding "higlier harmonic flows." The 
recipe: (a) project all 2D angular correlations onto ID azimuth, (b) perform 
a Fourier series fit to all data (which must describe any distribution on 
periodic azimuth), (c) interpret all Fourier series terms as "flows." In effect, 
the SS 2D jet peak is "fragmented" (via Fourier series) to become flows. 




Fig. 5. First: Correlations for 10-20% central Au-Au, Second: Projection with 
Fourier coefficients, Third: Comparing predicted multipoles with LHC data. 

Figure [5] (flrst, second) show 2D angular correlations from 10-20% central 
Au-Au (2D histogram) projected onto ID azimuth (solid points) and fitted 
with a four-term Fourier series (thin solid curves). The m = 3 (sextupole) 
Fourier term is interpreted as "triangular fiow" [13] . The sextupole structure 
(aka Mach Cones [7j) is derived from the jet-related SS 2D peak. 

Figure [S] (third) shows calculated Vm{SS} trends for SS 2D peak data 
from Ref. [5J (broken curves) and the nonjet quadrupole trend V2{2D} from 
Ref. [9] (solid curve). The predicted W2{2} (upper dotted curves) accurately 
describes published data from Ref. [Tlj. The points show a comparison 
of predictions with LHC data from Ref. [2]. The agreement is excellent 
modulo an overall factor 1.3 increase [8]. From the data we conclude that 
jet and nonjet-quadrupole trends are similar at RHIC and LHC energies. 
Consequences of adding a sextupole term to 2D model fits or applying such 
analysis below the Au-Au sharp transition are further described in Ref. [15] 

6. pQCD folding integral and fragment distributions 

We now demonstrate that the structure underlying claims of "higher 
harmonic flows" is described quantitatively by pQCD. To describe fragment 
yields in nuclear collisions requires a fragmentation function (FF) ensem- 
ble over the parton spectrum and the pQCD parton spectrum itself. The 
FF ensemble on dijet energy Q is accurately described (to statistical lim- 
its) by a beta distribution with rescaling of the total rapidity ^6j. FFs are 
represented by D{x,Q'^) D{y,ymax) = '^nchiymax)l3[u;p{ymax), q{ymax)]- 
Model parameters (p, q) are nearly constant over a large energy range. The 
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dijet fragment multiplicity ^nchiUmax) is simply determined by energy con- 
servation and scales with jet energy approximately as y^az = ln^((5/m,r)- 




Q = \'s(GeV) p, (GeV/c) 



Fig. 6. First: Measured fragmentation functions (FFs), Second: FFs rescaled 
and beta distribution, Third: FF ensemble vs energy Q, Fourth: pQCD parton 
spectra. 

Figure [6] (first, second) show measured FFs on rapidity y = ln[{p + 
E)/mTj] for three dijet energies and the same data rescaled on normalized 
rapidity u ~ y/Vmax [16j. For jet energies relevant to nuclear collisions more 
than half the jet fragments fall below pt = 1 GeV/c (y w 2.7), falsifying the 
assumption that hydro dominates below 2 GeV/c. The third panel shows 
the parametrized ensemble of FFs over a large energy range up to LEP II. 

Figure [6] (fourth) shows a parametrized parton spectrum (solid curve) 
compared to jet spectrum data (points) and pQCD theory calculations 
(other curves). pQCD describes the parton spectrum down to Q = 6 GeV (3 
GeV jets) corresponding to perturbative length scale 0.03 fm. The parton 
spectrum folded with FFs as the measured non-perturbative aspect of QCD 
determine /raymeni distributions (FDs) down to zero hadron momentum [6]. 
pQCD FDs can then be compared with spectrum hard components H [1]. 

7. pQCD descriptions of spectrum hard components 

Figure [7] (first) shows a calculated FD (solid curve, [B]) compared to 
spectrum hard component H derived from 200 GeV non-single-diffractive 
p-p data yy. A single parameter in the pQCD spectrum, the lower bound 
on the spectrum near 3 GeV, has been adjusted to accommodate the data. 
The data description is excellent. Note the mode of the FD near 1 GeV/c. 

In p-p collisions the two-component spectrum model is {1 / ns)dnch/ Vtdyt = 
Soivt) + {n'h/ns)HQ{yt) In A-A collisions the model is generalized to 
{2/Npart)dnch/ytdyt = SNNivt) + i^H{yt,i'), with nh/ris = 2Nun/Npart 
as analogous "centrality" parameters |4j. We expect hard component H{yt, u) 
to be modified in more-central A-A collisions due to modification of FFs. 
FF shape modification is modeled by changing parameter q in f3{u;p,q) by 
Aq. Fragment number 2nch is rescaled so as to conserve the parton energy. 
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y, 4p = lii(P|,/p) y = ln[(E + p)/m,] ^ y,(lt) 



Fig. 7. First: pQCD calculated fragment distribution (FD) vs p-p data, Second: 
Modified FFs on ^p, Third: Modified FFs on y, Fourth: pQCD FD vs Au-Au data 



Figure [7] (second, third) illustrates FFs for two dijet energies without 
(vacuum) and with (medium) modification. Figure [7] (fourth) shows the 
measured spectrum hard component for pions from 0-12% central 200 GeV 
Au-Au collisions (points, [4J) and the calculated FD with modified FFs 
(solid curve, [6j). Parameter Aq 1) has been adjusted to describe the 
five- fold suppression at yt = 5 {pt ~ 10 GeV/c). All else remains the same 
as for the first panel. The description of Au-Au data is remarkable. The 
dotted curve labeled GLS is the FD prediction for no FF modification. 



8. Resolved minijets and hadron yields vs Au-Au centrality 




Fig. 8. First: SS 2D peak volume, Second: pQCD calculated dijet frequency Third: 
Fragment yield vs Au-Au centrality Fourth: Total hadron yield vs centrality. 

Figure [8] demonstrates the correspondence between jet-related angular 
correlations and hard-component hadron yields with Au-Au centrality. The 
first panel shows the SS 2D peak volume within acceptance Arj = 2 (solid 
curve). The second panel shows the pQCD predicted dijet frequency within 
Arj. By combining the two trends the mean jet fragment multiplicity can 
be inferred. Recombining the fragment multiplicity with the dijet frequency 
predicts the hard component yields (solid curve) in the third panel and (com- 
bined with fixed soft component Snn) the total hadron yields in the fourth 
panel. The prediction (solid curve, [T7]) is compared with the measured 
hadron yields from Ref. [Ij (points). Again the agreement is remarkable. 
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9. Summary 

This presentation reviews quantitative relations among (a) jet-related 
2D angular correlations, (b) jet-related spectrum hard components and (c) 
pQCD-calculated jet frequencies and fragment distributions that strongly 
support a jet interpretation of the SS 2D peak for all Au-Au centralities. A 
pQCD description of hadron production at RHIC agrees with spectrum data 
within their uncertainties and implies that one third of final-state hadrons 
in central Au-Au collisions are contained within resolved jets. The more- 
peripheral 50% of the Au-Au total cross section corresponds to transparent 
collisions where dijet characteristics are just as for p-p collisions. Over 
the same transparency interval the NJ azimuth quadrupole measured by 
a statistically equivalent quantity increases to 2/3 of its maximum value, 
making a conventional hydro interpretation very unlikely. Proponents of 
"higher harmonic flows" must confront the many results that falsify a flow 
hypothesis and support a jet mechanism for the dominant correlation and 
spectrum structure. Cherry picking data features and favoring analysis 
methods and "theories" that seem to support a preferred (flow) hypothesis 
while disregarding contrasting results from pQCD is questionable practice. 
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